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Abstract 

We investigate the possibility of optical current control through single molecules 
which are weakly coupled to leads. A master equation approach for the trans- 
port through a molecule is combined with a Floquet theory for the time-dependent 
molecule. This yields an efficient numerical approach to the evaluation of the cur- 
rent through time-dependent nano-structures in the presence of a finite external 
voltage. We propose tunable optical current switching in two- and three-terminal 
molecular electronic devices driven by properly adjusted laser fields, i.e. a novel 
class of molecular transistors. 

Key words: molecular electronics, quantum control 
PACS: 85.65.-Fh, 33.80.-b, 73.63.-b, 05.60.Gg 



1 Introduction 



Spurred by the ongoing experimental progress in the field of molecular elec- 
tronics [1-5], the theoretical interest in transport properties of molecules has 
revived [6]. Tight-binding models for the wire have been used to compute 
current- volt age characteristics, within a scattering approach [7] and from elec- 
tron transfer theory [6] . Both approaches bear the same essentials [8] . For high 
temperatures, the wire electrons loose their quantum coherence and the trans- 
port is dominated by incoherent hopping between neighbouring sites [9]. Re- 
cently, the current- voltage characteristics has been obtained from a quantum- 
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chemical ab initio description of the molecule [10]. The results were in good 
agreement with recent experiments [5]. 

Typical electronic excitation energies in molecules are in the range up to an eV 
and, thus, correspond to light quanta from the optical and the infrared spec- 
tral regime where most of today's lasers work. It is therefore natural to use 
such coherent light sources to excite molecules and to study their influence 
on the transport properties aiming to find a way of manipulating currents. 
One particularly prominent example of quantum control is the so-called co- 
herent destruction of tunnelling (CDT), i.e. the suppression of the tunneUing 
dynamics in an isolated bistable potential by the purely coherent influence 
of an oscillating bias [11-15]. The crucial point there is that the long-time 
dynamics in a periodically driven quantum system is no longer dominated 
by the energies, but rather by the so-called quasienergies [15-17]. The latter 
may be degenerate for properly chosen driving parameters yielding a divergent 
time-scale. Inspired by these results, we address in this Letter the question of 
controlling by use of properly tailored laser fields the transport through time- 
dependent open systems, i.e. systems that allow for a particle exchange with 
external leads. 

For the computation of electrical currents through wires exposed to strong 
laser fields, we put forward Floquet approach [18]. The central idea of this 
method lies in a non-perturbative solution of the Schrodinger equation of 
the isolated time-dependent wire plus laser field, while the wire-lead coupling 
is treated pcrturbativcly. The resulting density operator equation is decom- 
posed into a time-dependent Floquet basis permitting a numerically efficient 
treatment. We generalise here this method to the analysis of networks with 
an arbitrary number of contacts to leads. Subsequently we apply the so ob- 
tained formalism to the investigation of optical current switching in two- and 
three-terminal devices as prototypical examples of a new class of molecular 
transistors. 



2 Model for wire and leads 

We embark by specifying the model Hamiltonian of the entire system as 
sketched in Fig. 1. It consists of the molecule in the laser field, ideal leads, 
and the molecule-leads coupling Hamiltonian, 

H{t) — -^^molecule(0 ~l~ -^leads ~l~ -^molecule-leads ■ (l) 

The irradiated molecule is modelled by a tight-binding description taking into 
account molecular orbitals \n), which are relevant for the transport. Disre- 
garding the electron-electron interaction, the most general form of the Hamil- 
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Fig. 1. Molecular circuit consisting oi N = 6 sites of which the sites 1, . . . ,L are 
coupled to L = 4 leads. 

tonian reads 

-f^molecule(^) = XI ^nn'{t) (^rf-n' ■, (2) 

where the fermionic operators c„ and c]j destroy and create, respectively, an 
electron in the molecular orbital \n). The sums extend over all tight-binding 
orbitals. The T-periodic time-dependence of the single-particle Hamiltonian 
Hnn'{t) = Hnn'it + T), reflects the influence of the laser field with frequency 
Q, — 2t:/T. The L ideal leads are described by grand-canonical ensembles of 
electrons at temperature T and electro-chemical potential I — 1,...,L. 
Thus, the lead Hamiltonian reads i^ieads = Z^gf cJ^Cg^, where Cg^ destroys 
an electron in state q in lead L All expectation values of lead operators can 
be traced back to {c\g^Cqie,i) = Sqq'Su'f{eq£ - m), where /(e) = (1 + e^/'^BTyi 
denotes the Fermi function. The model is completed by the molecule-leads 
tunnelling Hamiltonian 

-^^molecule-leads 

that connects each lead directly to one of the suitably labelled molecular 
orbitals. Since we are not interested here in the effects that arise from the 
microscopic details of the molecule-lead coupling, we restrict our analysis in 
the following to energy-independent couplings, i.e. = {27: /h) Y^q [Vqel"^ 5{e — 
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3 Perturbative description and Floquet ansatz 



Let us assume that the dynamics of the driven wire is dominated by the time- 
dependent wire Hamiltonian so that the couphng to the leads can be taken into 
account as a perturbation. This allows to derive by use of standard methods 
the approximate equation of motion for the total density operator Q{t) , 

Qit) = - ^[i^molecule(^) + ^^leads, Q{t)] (4) 

oo 

If — 
J~2 J d-T [i^molecule— leads ) [-^^molecule— leads '^i^^i ^(0]]- 



We have omitted a transient term that depends purely on the initial prepa- 
ration. The tilde denotes operators in the interaction picture with respect to 
the molecule plus the lead Hamiltonian, 0{t,t') = UQ{t,t')0{t)Uo{t,t'), where 
Uo{t,t') is the time-evolution operator without the coupling. The first term 
describes the coherent dynamics of the electrons on the wire, while the second 
term represents incoherent hopping of electrons between the leads and the 
wire. 



The net (incoming minus outgoing) electrical current that flows from lead i 
into the molecule is then given by the rate at which the electron number in 
the corresponding lead decreases multiplied by the electron charge — e, 

m^e^^iN,) . (5) 

Note that this expectation value is time-dependent through the non-equilibrium 
density operator g(t). To evaluate the right-hand side of Eq. (5), we employ 
Eq. (4) to derive after some algebra the result 

Re / 6-^/^(6 - f^e){[ce, 4{t - r, t)U) 

(6) 

- I ■ 

This expression still contains the yet unknown expectation values of solely 
those wire operators with a direct connection to lead £. It depends in par- 
ticular on the Heisenberg operators cj(t — r, t) and thus implicitly on the 
dynamics of the driven wire. Let us therefore focus on the single-particle 
dynamics of the periodically time-dependent wire Hamiltonian Hnn'{t)- An 
established procedure for the solution of the corresponding Schrodingcr equa- 
tion is to employ a Floquet ansatz amounting to a non-perturbative treatment 
of the external driving. There one uses the fact that a complete set of solu- 
tions is of the form \^>a{t)) = exp(— ieat//i)|$a(t)). The so-called quasiener- 
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gies Sa take over the role of the energy eigenvalues in static systems and 
govern the long-time dynamics. The Floquet modes obey the time- 

periodicity of the driving field which allows to express them as a Fourier 
series, \^a{t)) = Z^^L-oo ^^P(~'^^^^)l^a,fc)- They can be obtained from the 
eigenvalue equation [16,17,15] 

\n)H^n>{t){n'\-ih^)\^^{t)) = e^\^^{t)) . (7) 

n,n' ^'■^ 

Moreover, the Floquet modes define the complete set of operators 

Ca{t) = J2{^a{t)\n)Cn, (8) 
n 

whose time-evolution assumes the convenient form Ca{t—T, t) — exp{ieaT/h)ca{t). 
The orthogonality of the Floquet states at equal times [16,17,15] yields the 
back-transformation c„ = J2a{M^ai't)) ^-iid thus results in the required spec- 
tral decomposition. Using (8) and performing the energy and the r-integration 
in Eq. (6), we obtain for the time- averaged current the main result 



E 



^ ak 



{^^,k\e){e\^^,k)f(eo. + khn-iie) 

n (9) 

0k' 

Here, we have introduced the expectation values Rapit) = {Ca{t)ci3{t)) , which 
assume in the long-time limit the time-periodicity of the driven system and 
thus can be decomposed into the Fourier series Rapit) = J2k ^w{~'ikfit) Rai3,k ■ 
It is straightforward to derive for the Ra/3,k from the density operator equa- 
tion (4) the following set of inhomogeneous linear equations 

ji^a - e/3 + knn)Rap,k (10) 

= qE^^I ^{^l3,k"+k'\^){^\^P',k"+k)Ral3',k' 
^ ek' ^ I3'k" 

+ E {^a',k"+k'\^){^\^a,k"+k)Ra'l3,k' 
a'k" 

- {^p,k'-k\i){i\^a,k')f{ea + k'n - fie) 

- {<^>p,k'\i)ma,k'+k)fief, + k'n - fl,) 

which will be solved numerically. We have found that even in the case of strong 
driving where the Floquet states comprise many sidebands, a few Fourier 
coefficients Rai3,k are in fact sufficient to obtain numerical convergence. This 
justifies a posteriori the use of the Floquet states as a basis set. 

To conclude the technical part of this work, we note that our approach goes 
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beyond a linear response treatment of the driving and additionally does not 
use a so-called rotating wave approximation (RWA) [19,20], which neglects 
the oscillatory contributions to Ra(}{t) by the ansatz Rap,k = Pa^a,i3^k,o- 

In 

fact, we found that a RWA solution dehvers inaccurate results in the vicinity 
of quasienergy degeneracies. 



4 Optical current gate 

As a first setup that may be suitable as a current control device, we investigate 

the transport through a two-level system, i.e. a wire that consists of = 2 
sites — one of them is coupled to the left lead and the other to the right lead. 
Then the time-dependent wire Hamiltonian reads in the basis of the molecular 
orbitals 



where A denotes the tunnel matrix element between the two sites and E'l 
and £'r are the corresponding on-site energies. The laser field contributes to 
the Hamiltonian (11) a time-dependent bias with amplitude A — —e£d, i.e. 
charge times electrical field strength times the site-to-sitc distance. Note that 
the electrical field may be drastically enhanced due to the presence of the 
metallic tips [21]. The effective coupling to each lead is assumed to be equal, 
= r, and an external voltage V is taken into account by a difference in the 
electro-chemical potentials, /^l — A*r = —eV. 

We use in all numerical calculations the tunnel matrix element A as the energy 
unit and assume that the effective couplings to the leads are by one order of 
magnitude smaller, hV — O.IA. This corresponds to a large contact resistance 
and ensures the applicability of a pcrturbational approach. A realistic value 
is A = O.leV, resulting in a current unit eV = 0.256 /xA. For a sitc-to-site 
distance of lOA and a laser frequency il = 10 A/^, the driving amplitudes 
considered below correspond to an electrical field amplitude of lO^V/cm at 
1 /im wavelength. 

The time-averaged current I = Ij^ = — Jr through the molecule in a case 
where both on-site energies are equal is depicted in Fig. 2. As a striking fea- 
ture, we find that at certain values of the driving amplitude, the current col- 
lapses to less than 1% of its maximal value reached in the absence of the 
driving. Closer inspection (not shown) reveals that the width of this depres- 
sion is proportional to the molcculc-lead coupling F. Comparison with the 
quasienergy spectrum in the lower panel demonstrates that the current break 
downs occur at quasienergy crossings. This relates the present phenomenon to 




(11) 
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Fig. 2. Average current and quasienergy spectrum versus driving amplitude for 
a wire which consists of two sites between two electrodes (cf. inset) for unbiased 
(E'r = £'l = 0, solid lines) and biased (£^r = — £'l = 0.1 A, dashed lines) on-site 
energies. The leads' chemical potentials are /xr = — /xl = lOA; the other parameters 
read hn = lOA, ksT = 0.25A, HT = O.IA. 

the CDT, i.e. the standstill of the tunnel dynamics in a driven bistable poten- 
tial at quasienergy crossings [11]. For the isolated two- level system (11) with 
A <S A, the CDT condition reads Jo(2A/^f2) = [12], i.e. the suppression 
of the tunnelling dynamics is related to the zeros of the Bessel function Jq. 
As our analysis reveals, the same condition results in a suppression of the 
transport through the open system. 

An external voltage may be of peculiar influence to the on-site energies of a 
molecular wire [22] and cause an effective bias -El 7^ -E'r in originally symmetric 
molecules. Thus, a crucial question is whether the above current suppression 
is stable against such a modification. The broken lines in Fig. 2 demonstrate 
that this is indeed the case. Although the quasienergies now form an avoided 
crossing, the current breakdowns do survive; they are even more pronounced, 
but slightly shifted towards larger driving amplitudes. This robustness of CDT 
based current control combined with the huge on/off ratio suggests the pre- 
sented setup as a promising alternative to structural chemistry-based switch- 
ing devices [23,24]. 
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Fig. 3. Average currents through contacts Ci (sohd) and C2 (broken) as a function 
of the polarisation angle (f) for the three-terminal device depicted in the inset. The 
chemical potentials are /xe = — A*Ci = — MC2 = 50A; the on-site energies = 0. 
The driving field is specified by the strength A = 25A and the angular frequency 
O = lOA/h; the effective coupling is HT = O.IA and the temperature ksT = 0.25A. 
The maximal value of the current ratio /ci//c2 ~ 100 is assumed &t (j) = 60°. 

5 Optical current router 

An experimentally more ambitious configuration consists in a planar molecule 
with = 4 sites, three of which are coupled to a central site and are directly 
connected to leads (cf. inset of Fig. 3). We borrow from electrical engineering 
the designation E, Ci, and C2. Here, an external voltage is always applied such 
that Ci and C2 have equal electro-chemical potential, i.e. /iCi — A*C2 7^ A*e- In 
a perfectly symmetric molecule, where all on-site energies equal each other, 
reflection symmetry at the horizontal axis ensures that any current which 
enters at E, is equally distributed among Ci,2, thus Ic^ = = — /e/2. Since 
this structure in Fig. 3 is essentially two-dimensional, we have to take also 
the polarisation of the laser field into account. We assume it to be linear 
with an polarisation angle as sketched in the inset of Fig. 3. The effective 
driving amplitudes of the orbitals which are attached to the leads acquire now 
a geometric factor which is only the same for both orbitals Ci and C2 when 
(p — 0. For any other polarisation angle, the mentioned symmetry is broken 
and the outgoing currents may be different from each other. The difference 
may be huge, as depicted in Fig. 3. Their ratio varies from unity for = up 
to the order of 100 for = 60°. Thus, adapting the polarisation angle enables 
one to route the current towards the one or the other drain. 

Alternatively, one can keep the polarisation angle at ^ = and break the 
reflection symmetry by using an intrinsically asymmetric molecule, as sketched 
in the inset of Fig. 4. This allows to control sensitively the ratio of the outgoing 
currents by the strength A of the external field, cf. Fig. 4. The switching range 
comprises up to four orders of magnitude with an exponential sensitivity. 
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Fig. 4. Ratio of the outgoing average currents versus driving strength A for the 
three-terminal device at a polarisation angle ^ = 0. The filled circle in the inset 
depicts a site with an on-site energy Eq^ that differs from the others. All other 
on-site energies and parameters as in Fig. 3. 



6 Concluding remarks 



We have presented a method for the efficient numerical computation of cur- 
rents through periodically time- dependent networks with two or more contacts 
to external leads. The application to two types of setups substantiated that 
external fields bear a wealth of possibilities for the manipulation of electrical 
currents: in a molecular wire the current can be suppressed by proper time- 
dependent fields. In a three terminal device, it is possible to route by tailored 
optical fields the current that enters from a source towards the one or the 
other drain. 

The authors hope that their proposals will motivate experimentalists to accept 
the challenge of implementing the proposed molecular wire schemes in the lab- 
oratory. The two-terminal current gate can possibly be realized using break 
junctions exposed to a laser field. Alternatively, one could use a self-assembled, 
laser-irradiated maze-like layer of sparsely distributed conducting molecules on 
a conducting surface. Then by positioning a scanning tunnelling microscope 
tip directly over one such molecule, it should be possible to measure the fea- 
tures of the predicted gating behaviour. Experimentally more ambitious is the 
realization of the arrangement in Fig. 3 with a planar molecule contacted to 
three electrodes. Here again, laser-irradiated self-assemblies of molecules such 
as carbon nanotube complexes or of biomolecules hke metalised DNA [25], or 
the use of cationic hpid-DNA complexes [26] as DNA-nanocables, with the 
centre-molecule covalently bound to such planar structures, might make the 
experiment feasible. 

A completely different realization of our findings should be possible in semi- 
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conductor hcterostructures. There, instead of a molecule, coherently coupled 
quantum dots [27] form the central system. Furthermore, owing to the lower 
level spacings, the suitable frequency of the coherent radiation source is then 
in the microwave spectral range. 
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